X 



Dissipative structures in a 
nonlinear dynamo 



Andrew D. Gilbert,^ Yannick Ponty^ 
and Vladislav Zheligovsky *^ 

^^ '^Mathematics Research Institute, University of Exeter, U.K. 

(^ Observatoire de la Cote d'Azur, Nice, France. 

^International Institute of Earthquake Prediction Theory 
^ and Mathematical Geophysics, Moscow, Russia. 

1—5 

o 

en 

^ Abstract 
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\^ This paper considers magnetic field generation by a fluid flow in a system referred 

r^ to as the Archontis dynamo: a steady nonlinear magnetohydro dynamic (MHD) state 

'rl^ is driven by a prescribed body force. The field and flow become almost equal and 

, ", dissipation is concentrated in cigar-like structures centred on straight-line sepa- 

ratrices. Numerical scaling laws for energy and dissipation are given that extend 

\jj previous calculations to smaller diffusivities. The symmetries of the dynamo are 

qn^ set out, together with their implications for the structure of fleld and flow along 

•/^ the separatrices. The scaling of the cigar-like dissipative regions, as the square root 

Sp^ of the diffusivities, is explained by approximations near the separatrices. Rigorous 

• results on the existence and smoothness of solutions to the steady, forced MHD 

in ^. . 

r^ equations are given. 
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1 Introduction 



Much is known about fast dynamo action: the rapid growth of magnetic fields 
at higli magnetic Reynolds number in fluid flows with chaotic streamlines, 
but the mechanisms for the dynamical saturation of such fields remain poorly 
understood. In many cases when the growing field equilibrates by modifying 
the fluid motion, the effect is to switch off the chaotic stretching in the flow, 
as measured for example by a reduction in the finite-time Liapunov exponents 



e.g., Cattaneo, Hughes & Kim, 1996 Zienicke, Politano & Pouquet 1998). 



What is left is a fluid threaded by a magnetic field which resists stretching 
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and so suppresses overturning fluid motions, but supports elastic wave-like mo- 
tions, essentially Alfven waves with coupled fleld and flow (e.g., Courvoisier. 



Hughes & Proctor, 2010). The flnal state of many simulations shows appar- 



ently chaotic behaviour in space and time, suggestive of an attractor of moder- 
ate or high dimension, although because of the three-dimensionality of MHD 
systems little can be done to explore its properties, for example the fractal 
dimension or spectrum of Liapunov exponents. 



Although this appears to be the outcome of many simulations, as far as they 
can be run, there are some intriguing examples where a further phase of evo- 
lution takes place: the magnetic fleld and flow align, depleting the nonlinear 
terms, and both flelds evolve to a steady (or very slowly evolving) state. The 
key point is that in unforced, ideal magnetohydrodynamics (see equations 
(2.1-2.3) below with u = r] = and / = 0) any state with u = ±b is an 



exact steady solution. The remarkable fact that simulations of forced, non- 
ideal MHD turbulence could evolve to something very close to such a state 
was first observed by Archontis ( 2000[ ) in his thesis, and published in Dorch 
& Archontis (2004) (hence referred to as DA), and Archontis, Dorch & Nord- 



lund (2007). These simulations use a compressible code with a Kolmogorov 



forcing function, (2.4) below, first used as the form of a fiow for simulations of 



fast, kinematic dynamo action by Galloway & Proctor (1992). Subsequently 



Cameron & Galloway (2006a) undertook incompressible simulations of the 



same system as Archontis, and pushed up the fiuid and magnetic Reynolds 
numbers; our work is linked closely to this paper, which we refer to as CG in 
what follows. 



What these authors found was that, starting with a forced fiuid fiow and a 
seed magnetic field, the growing magnetic field initially equilibrates in rough 
equipartition with the velocity field, in a messy, chaotic time-dependent state. 
However during this state, there is a slow but persistent exponential growth in 
the average alignment of the u and b vectors, as measured by the cross-helicity. 
This process of alignment continues until there takes place a sudden increase in 
the fiuid and magnetic energies, and both fields tend to a steady state of almost 
perfect alignment, discrepancies being controlled by the weak dissipation and 
the forcing. In fact since any solution u = ±b is a neutrally stable solution of 



the ideal problem (Friedlander & Vishik, 1995), the solution that is selected 



must depend delicately on balances involving these sub dominant diffusive and 
forcing effects. We note that some alignment of field and fiow has been noted 
in many other MHD fiows, for example see IDobrowolny, Mangeney & Veltri 



1980 



Pouquet, Meneguzzi &; Frisch (1986), Mason, Cattaneo & Boldyrev 



(2006) and references therein, but of a less dramatic nature. 



This observation of dynamo saturation in a steady state with such a high 
degree of alignment was a new phenomenon: CG refer to the saturated state 
as the 'Archontis dynamo', though we prefer the term 'Archontis saturation 



mechanisin' . CG observed this aligned state as a solution branch over a wide 
range of magnetic and fluid Reynolds numbers (taking the magnetic Prandtl 
number to be unity in much of their work). Further developments include 
the development of bursts of rapid time dependence after some time in the 



steady state, in the study Archontis, Dorch & Nordlund (2007). However this 



appears only to occur in the compressible case, as it has not been seen by 



CG nor in our simulations; we will therefore not discuss this further. Cameron 



& Galloway (2006b) also find slow time-dependent evolution of the saturated 
state for the Kolmogorov forcing with magnetic Prandtl number Pr = v/rj 
not equal to unity, and for more general spatially periodic steady forcings. In 
all cases though, the field and flow settle into a state of very close alignment, 
even if they then evolve on a slow time scale. 



The focus of the present paper is to understand more about the structure 
of the steady saturated state for the Kolmogorov forcing and unit magnetic 
Prandtl number Pr, with a particular focus on the regions where dissipation 
occurs and on rigorous results on existence and smoothness. DA and CG find 
a complex geometrical picture for the field and flow and identify these regions 
of high dissipation: they are localised along straight-line separatrices that join 
a family of stagnation points; similar structures are found in the 1:1:1 ABC 



flow (Dombre et al, 1986). These are found to have a width scaling as ^/e 
where £ is a dimensionless measure of the diffusivity, and one of our aims is 
to understand this power law. 



We set up the governing equations in ^ and extend the solution branch to 
yet smaller values of the diffusivity e by means of large scale simulations in 
^ In ^ we then classify the symmetries of the Kolmogorov forcing, which 
are preserved by the nonlinear, saturated field and flow. These symmetries are 
the reason for the presence of the non-generic straight line separatrices that 
join stagnation points in the flow and field, and they constrain the local flow: 
it is in these regions that dissipation is strongest. We plot the local structure 
of fields along the separatrix from (0, 0, 0) to (tt, tt, tt) in ^ We determine the 
effects of diffusion by setting up PDEs for the advection of field as it enters 
the dissipative regions in ^ and use these to justify the order ^/e scaling for 
the cigar widths found in CG. We then proceed with a formal mathematical 
investigation of the existence of steady-state solutions to the MHD problem 
at hand and bounds for them in various function spaces in §^|^-[9j The reader 
should note that these sections use functional analysis and so have a different 
flavour from the earlier ones. Finally ^|T0] offers concluding discussion. 



2 Governing equations 



We begin with tlie dimensional equations for incompressible MHD, in the form 

dtU + uVu = bVb-Vp + uV'^u + f, (2.1) 

dtb + uVb = bVu + T]V^b, (2.2) 

V-u = V-b = 0, (2.3) 

where z/ and t] are the kinematic viscosity and magnetic diffusivity. We take 
/ to be a steady force of magnitude J-" acting on a length scale C. We will 
consider the Kolmogorov forcing f = J^ f* [r / C) , whose dimensionless form is 
given by 

f*{r) = {sin z,smx,smy). (2.4) 

In non-dimensionalising we have only the parameters {£, J-", z/, r^}, together 



with the form (2.4) of the forcing function. From these we can define a mag- 
netic Prandtl number and a Grashof number as in similar forced flow problems 
(see, e.g., 



Childress, Kerswell & Gilbert, 2001) by 



Pr = u/t], Gr = J^C^u' = e-\ (2.5) 

We have as diagnostics the Reynolds number and magnetic Reynolds number 
given by 

Re = C\\u\\/u, Rm = C\\u\\/ri, (2.6) 

where \\u\\ is a measure of the fluid velocity at a given time, for example the 
L^ norm, taken as the root-mean-square value, averaged over the periodicity 
box. We rescale as 

u = Uu\ b = Ub*, t = {C/U)t*, r = Cr\ f = J^f\ p = U^p*, 

(2.7) 
with the choice of velocity scale 

U = TC?lv. (2.8) 

This yields the non-dimensional formulation, dropping the stars, as 

9tw + w- Vu = b- Vb- Vp + eV^u + e/, (2.9) 

9*6 + w ■ Vb = b • Vw + ePr" V^b, (2. 10) 

V-u = V-b = 0, (2.11) 



with / given in (2.4) and the only parameters specified are {e,Pr}. The cor- 



responding Reynolds and magnetic Reynolds numbers are 

Re = £"^11^11, Rm = £-^Pr||'u||. (2.12) 



We refer to e~^ as the Grashof number Gr and will be interested in the in- 
viscid limit e — )■ 0. The Reynolds number and magnetic Reynolds number 
are diagnostics depending on the flow regime realised^ Indeed, they change 
greatly during the saturation process, when the fields align and 11^11, \\b\\ in- 
crease significantly. As in CG, the governing equations may be written in a 
more symmetrical form in terms of Elsasser variables 



u±b, (2.13) 



which gives, for Pr = 1 



dt\+ + \--V\+ = -Vp + e\/^\+ + Ef, (2.14) 

9iA_ + A+ ■ VA_ = -Vj9 + eV^A^ + e/, (2.15) 

V-A+ = V-A_=0. (2.16) 



3 Numerical results 



We undertook a number of runs to investigate the structure of the steady, 
equilibrated Archontis dynamo for Pr = 1 and values of e down to 10~^ in the 
(27r)^ periodic domain T^. The steady solutions were found by following the 
solution branch: that is taking the output from a run with a given value of e 
and using it as the initial condition for a run with a reduced value of e. This 
establishes the Archontis dynamo as a robust local attractor, in the range 
of e used, in agreement with DA and CG. Whether it is a global attractor 
over some or all sufficiently small values of e remains unknown, and extremely 
difficult to address in view of the long transients that may occur. Our runs 
were undertaken with a pseudo-spectral code using A^^ modes with A^ = 128 
for e = 0.02 and 0.01, A^ = 256 for e = lO'^, and A^ = 512 for e = 10^1 
There were other, less well resolved runs with A^ = 128 for e = 10"'^ and 
A^ = 256 for e = 10"'*, which we refer to below as our 'testing simulations'. 
For comparison, CG go down to £ = 1.25 x 10~^ in their study, with resolution 
128^. Our results thus extend theirs by a little over a decade, and in this section 
we present measures of the magnetic field and flow in the equilibrated state. 

Numerical values are given in table 1 and plotted in flgure[lj Panel [I]^a) shows 
the kinetic and magnetic energies in the equilibrated state, given by 

Ek= f l\u\^dV, Em= f l\b\^dV. (3.1) 



^ Our formulation is equivalent to DA/CG, but our terminology is a little different. 
For example CG use the parameters vcq = e and r/cc = ePr~ , which they refer 
to as the inverse Reynolds and magnetic Reynolds numbers respectively. 
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Fig. 1. Numerical results plotted against logj^g^^^- Plotted are (a) kinetic energy 
Ek (triangles) and magnetic energy Em (squares), (b) enstrophy Qk (triangles) 
and squared current Qm (squares), (c) normalised cross helicity Hx/C^EkEm), (d) 
energy log^o -S- of A_ (dotted line gives e^ dependence). 
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N 


Ek 


Em 


^K 


^M 


Hx 


E^ 


0.02 


64 


0.1797 


0.1745 


0.1849 


0.1791 


0.3532 


8.685 X 10"^ 


0.01 


64 


0.1781 


0.1765 


0.1825 


0.1816 


0.3543 


2.313 X 10-4 


0.001 


256 


0.1717 


0.1717 


0.1786 


0.1789 


0.3435 


3.04 X 10-*^ 


0.0001 


512 


0.1722 


0.1722 


0.1787 


0.1787 


0.3443 


3.55 X 10-s 



Table 1 
Numerical results. 

These show an initial decrease with e (as in CG) but then a shght increase 
from e = 10~^ to e = 10"^: this is quite small bearing in mind the scale on the 
vertical axis, but appears to be real as it is borne out in our test simulations. 
In all these runs Ek > Em though this is not apparent from the numbers 
in table 1 nor in panel [Ila). Panel fiyb) shows the enstrophy and integrated 
squared current, defined by 



Q 



K 



I \\vxu\^dv, nM= I \\vxh\^dv. 



(3.2) 



The total dissipation is given by 2eVLK + '2^s^m and this tends to zero as 0{e), 
as does the input of mechanical energy. Panel fltc) shows the cross helicity 



H 



X 



u-bdV 



(3.3) 



in normalised form, which rapidly tends to its theoretical upper bound of 
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Fig. 2. Numerical results for the case e = 10~^ with N = 256. Plotted are (a) IOEk 
(upper) and WEm (lower), (b) lOilj^ (lower) and IOQm (upper), against f/10. (c,d) 
are the same but plotted against 10^/(t — 1400). In each panel dotted lines are linear 
fits. 

unity, within the accuracy of our simulations, indicating the strong ahgnment 
of field as e — )■ 0. Finally panel [lid) shows the energy in the A_ Elsasser 
variable, where 



E, 



l\A±\'' dV = Ek + Em ± H^ 



X- 



(3.4) 



This shows a rapid decrease to zero as e — )■ consistent with the scaling 
E_ oc e"^ (dotted line) in agreement with the discussion in CG and below p] 

In panel [lib) it is notable that the two curves, for enstrophy and total current 
squared, cross between e = 0.01 and 0.001. The enstrophy Qk is a little 
smaller than Qm for e = 10~^ and in fact is also for 10~^ and in our test 
simulations, making us confident that this is a real effect. This opens up the 
question of how we measure these quantities, since the rate of evolution of the 
state becomes extremely slow for small e. Figure [2]^a,b) shows E^, Em, ^k 
and Qm as functions of time for the case e = 10~^ and A^ = 256: comparison 
with linear fits (dotted) shows clear curvature, as expected, but also highlights 
the slow evolution. This suggests neutral stability of the final state, and an 
expansion for any quantity in the form 



A = Ao + Ait-^ + A2r^ + 



(3.5) 



^ The (downwards) slope of a line fitting the data points is close to 1.9: the reason 
for the discrepancy is unclear: it could be numerical, or the e^ power law may only 
be achieved as e — )• which is possible as there is some downwards curvature present 
in the data points. 




Fig. 3. Cross sections showing |A_| in the (x,y)-plane for z = 7r/2. In (a) e = 10 ^, 
(b) 10~^ and (c) 10~^ and the colour scale shown runs from zero (bottom) to (a) 
0.031, (b) 0.0041 and (c) 0.00058 (top). 

Although asymptotically the origin of time does not matter, we found it helpful 
to choose an origin of time to (once per run) so as to obtain the best linear fit 
for quantities in the form 



A^v4o + Ai(t-to 



(3.6) 



We then use an estimate of the limiting value as Aq] for example see figure [2] 
(c,d). This was done for all the results in table 1. 

One of the aims of this paper is to focus on dissipative regions in the system: 
these occur along a series of straight line separatrices (DA/CG) and in figure 
[sl we show colour plots of |A_| for a range of diffusivities. Clearly seen in each 
case, but especially in (c) at the smallest e, are cross sections of spiralling field, 
centred on the separatrices, where small scales are generated with consequently 
enhanced diffusion. 



4 Symmetries 



We have seen that the dissipation tends to concentrate in cigar shaped regions, 
with one extending from (0,0,0) to the stagnation points at ±(7r,7r,7r). The 
reason these straight line separatrices are robust structures is linked to the 



study by Galloway & Proctor ( 1992 ) of the Kolmogorov flow 



symmetries of the forcing (2.4) and also applies to the kinematic dynamo 



■iiKoil^) = (sin 2;, sin X, sin y). 



(4.1) 



These symmetries turn out to be preserved by the solution (w, b) in the non- 
linear regime: there is no symmetry breaking. The forcing ( |2.4 ) is 27r-periodic 
in each coordinate: we only consider symmetries up to this periodicity (and 
that do not reverse time). Note first that any map T maps a vector field u 
according to 



{Tu){t) 



'T 



u{T 



-1„ 



drjdT 



-1 



r. 



(4.2) 



It is then easily checked that the forcing / in (2.4) is preserved by the fol- 
lowing 12 orientation-preserving symmetries, with det J = 1, which form the 
group y44 of even permutations of 4 objects, or the symmetry group of the 
tetrahedron, 



i{r) 


= ix,y,z), 


a\r) 


= i-x,7T-y,z + 7i), 


b\r) 


= ix + 7T,-y,7r-z), 


c'ir) 


= iT^ -x,y + 7T,-z), 


d{r) 


= iz,x,y), 


d\r) 


= iy,z,x), 


e(r) 


= i-z,TT -x,y + TT), 


2/ \ 

e (r) 


= in -y,z + Tr,-x), 


fir) 


= {z + n,-x,n -y), 


fir) 


= i-y,-K - z,x + tt), 


9ir) 


= i7T-Z,X + TT,-y), 


9» 


= iy + TT, -Z,TT - x). 



(4.3) 



These also form a subgroup of the group of 24 symmetries of the 1:1:1 ABC 



flow (Arnold & Korkina, 1983; Dombre et al. , 1986), and the above follows the 



notation in Gilbert (1992). The symmetries all commute with the inversion 
symmetry j{r) = (— x, —y, —z) and so the full symmetry group of the forcing 
/ is the direct product A4 x Z2. 



5 Flow and field on the separatrices 



The above symmetries constrain the behaviour of the magnetic field and flow 
on the separatrices. Take, for definiteness, the separatrix joining (0, 0, 0) to 
(tt, tt, vr) and call this the 'main separatrix' for brevity. Because of the sym- 



metries d and d in (4.3) there is a three-fold rotational symmetry about this 



separatrix, as seen in DA/CG, and any vector field on the separatrix can only 
point along the separatrix. We may introduce rotated Cartesian coordinates 



via 



/ \ 


/ 


X 


= 


w 


\ 



1/v^-l/v^ 
l/v^ 1/V3 1/V3 



4^ 



VJ 



(5.1) 



with ( along the separatrix. From there we may further define cylindrical polar 
coordinates (p, 6, Q, whose axis is ( along the separatrix with /i = pcosO and 
X = psinO. 

Our aim now is to investigate more of the behaviour of the flow near to the 
separatrix, in the saturated regime. However to fix ideas and establish a bench- 



mark, we consider first the Kolmogorov flow ■ukoI in (4.1). For this flow it can 



be shown that on the main separatrix motion is governed by 



C = V3sin(C/y3), p 



V 



(5.2) 




Fig. 4. Plot of components of fields against C for runs with e = 0.02, 0.01, 10^'^ and 
10~^ for (a) the field A+, reading down the curves (separated by adding 0, —0.1, 
—0.2, etc.), (b) the field V x A+, reading down the curves (c) the field e~-^/^V x A_, 
reading down the curves (separated by adding 0, —0.25, —0.5, etc.) (d,e,f) the field 
e-^A_ for e equal to (d) 0.01, (e) lO'^ and (f) 10"^. 



with solution 



C, = v3 (vr — cos ^ tanh t) . 



(5.3) 



Here ( tends to zero as t — )■ — oo and to a/Stt as t — > oo. Near to the separatrix, 
the radial coordinate p ^ 1 and the flow field may be expanded in powers of 
p. In view of the three-fold rotational symmetry, the flow u is axisymmetric 
about the main separatrix p = at leading order and streamlines are given by 



V(c)p + o(/), e = n{0 + o{p), c = 2s{0 + o{p'), 



(5.4) 



with 



2sKoi(C) = v^sin(C/y3), 2fiKoi(C) = y3cos(C/\/3). 



(5.5) 



Trajectories spiral in for C — and spiral out for ( ~ \/37r. On the separatrix 
itself u = (0, 0, 2s(C)) and V x w = (0, 0, 2Q{()), directed along the axis. 

Now in the nonlinear, equilibrated regime, the symmetries of the system are 
observed to be preserved and so the motion near and along the separatrix is 



given by (5.4) for some functions s{Q and fi(C). These functions characterise 



aspects of the nonlinear saturation on the separatrices and so of the spiral 
dissipative structures that form there, visible in figure [3j We can measure the 
equivalent functions for any field, and in our runs we find that the traces for 
2u, 2b and A+ are identical to graphical accuracy. In figure |4](a) we show the 
components of A+ along the separatrix (separated by constants). This figure 
in fact depicts two separatrices, the main separatrix from (0,0,0) to (7r,7r, tt) 
and the next one that continues to (27r, 27r,27r), with < C < 2\/3ti. The 
components of A+ show a sinusoidal form in keeping with the property of the 
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equilibrated fields noted by CG, namel}p] 

A+ ~ UKoh (5.6) 

There is only slight steepening at (vr, vr, vr) as e is reduced. Panel |4]^b) shows 
traces of the components of V x A_|_ with clear cosine form, in keeping with 



(5.6) and (5.5) but of somewhat enhanced amplitude, and with evidence of 



some finer scale structure near (tt, tt, it). These indicate that the approximation 



(5.6) is reasonable for the leading order fields on the separatrices. 



The picture is naturally more complicated for the A_ field, which tends to zero 
in the limit of small e. Panels E[d,e,f) plot the components of e~^A_ along 
the separatrix: there is clear evidence of finer scale oscillations emerging in the 
limit, but the nature of the limiting distribution is unclear. Panel |4Fc) shows 
the fields e~^^'^'V x A_ (separated by constants). These show the development 
of a cusp at (tt, tt, tt), the stagnation point where the two separatrices converge. 
In conclusion, the field A_ on the separatrix scales as 0{e), but its curl scales 
as 0(£^/^), giving a natural O^e^^"^) cigar width length scale, confirming results 
in CG and to be explored further below. 



6 Local behaviour and scaling in the cigars 

We now have some knowledge of the local structure of the fiow and field on the 



separatrices, in terms of both the general form it must take, namely (5.4), and 
the actual behaviour for small values of e in figure |4j The aim of the present 
section is to derive the dissipative lengthscale of y/e noted by CG. Of course 
we are not able to put together a solution that is complete: the dissipative, 
cigar-like regions process field that is drawn in, in a spiralling fashion, and 
then churn it out again. A complete picture would involve matching to the 
outer region, which is a highly three-dimensional problem, beyond what we 
can do; nonetheless a local picture gives some information. 



6.1 Uncurling the induction equation 



We start with the formulation in Elsasser variables (2.13-2.16) and for brevity 
set 

A = A+, £:A = A_, p ^ ep. (6.1) 



^ As noted by CG, although this is a good approximation, the error does not go to 
zero with e (e.g., |||A+ — ■ukoi|||oo/|||A+|||oo remaining at a level of about 15% for 



all runs) and so (5.6) should not be seen as an asymptotic statement. 
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We assume the key scaling of CG that A = 0(1), at least in the outer region, 
which means away from the stagnation points and the separatrices. Without 
approximation, the steady equations are 



A VA 

A- VA 



-Vp + V^A + Z, 
-Vp + eV^X + f. 



(6.2) 
(6.3) 



Note that a straightforward estimate of the width of a diffusive layer based on 
(6.3) would suggest an order e scaling from balancing A- VA ~ eV^A, but this 



is too small, as it does not take into account the different scales of variation of 
A along and across the characteristics of A, and the following, more delicate 
argument is needed. 



equation (2.10) 



Subtracting (|6.3|) from (|6.2|) gives an equation equivalent to the induction 

(6.4) 



= V X (A X A) + V^A - eV^A, 
which may be uncurled as 

Va = AxA-VxA + £:VxA, 



(6.5) 



where a{r) is a scalar field. Taking the divergence gives an elliptic equation 
for a, 

V^a = V ■ (A X A). (6.6) 

This development can be pursued further, to obtain a general closed but com- 
plicated system of scalar PDEs that link the field and flow to the external 



forcing, as in Zheligovsky (2009). However our present aims are more limited: 



we only need that (6.5) is equivalent to two equations. 



A-Va = -A-VxA + £A-VxA 



and 



A = cA + A^^A X (Va + V x A - eV x A), 
where c{r) is another scalar field which obeys 

A . Vc = -V ■ [A^^A X (Va + V X A - £V X A)], 



(6.7) 
(6.8) 

(6.9) 



from requiring that V ■ A = 0. Everything is exact up to this point but we 
note that this representation will generally break down at isolated points where 
A = 0. 

Now we approximate: first consider an 'outer' region, well away from the 
dissipative, cigar-like structures that lie on the separatrices joining stagnation 
points. We neglect diffusion in the outer region, the fields having a greater 
length-scale. The leading order outer problem is obtained by simply setting 
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e = in the above equations (6.7 6.9), leaving a pair of quasi-linear equations 



for a and c giving transport along characteristics of A, namely 

A-Va = -A-VxA, (6.10) 



A . Vc = -V ■ [A^^A X (Va + V X A)], 



(6.11) 



together with an equation that then reconstructs A, from (6.8), which we write 
as a sum of three terms, 



with 



Aa + Ay 



Ac = cA, A„ = A^^A X Va, Aa = A^^A x (V x A). 



(6.12) 
(6.13) 



Finally for this section, we note that in the outer region, Ac can be calculated 



explicitly in terms of A and a. Substitution of (6.13) into (6.3), where the 
diffusive term involving e is neglected, yields 



Vx 



(A ■ V) (cA + A-^A X (Va + V X A))] = V x /. 



By virtue of (6.11), this equation takes the form 



Vc X (A ■ V)A + cV X ((A ■ V)A) = V x F, 



(6.14) 



(6.15) 



where 



AV-[A-"Ax(Va + VxA)]-(A-V)[A-"Ax(Va + VxA)] + /. (6.16) 



Scalar multiplication of (6.15) by (A ■ V)A yields 

(V X F) • (A ■ VA) 



(6.17) 



[V X (A ■ VA)] • [A ■ VA] ■ 
Thus singularities of c can arise, where the helicity type term involving (A- V) A 



(i.e., the denominator in (6.17)) vanishes. 



6.2 Field in the outer region, near the main separatrix 



In the outer region, as the main separatrix is approached, it is observed that 
the field A is relatively smooth, as seen from the numerical simulations of CG, 



and also in view of the leading order Laplacian in (6.2), whereas A develops 



fine scales. Using the formulation in (5.4) we from now on define s{Q and ^{C) 
by 

A = i-s'iOp, Q{C)p, 2s{C)) + 0{p'), (6.18) 
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Fig. 5. Schematic figure showing the flow in the (/9, 9, Q) coordinates. 



in cylindrical polar coordinates {p,6X) defined in (5.1) and below. Here the 
functions 2s{Q/y/3 and 2^{Q\'^ defined for A are shown in figure |4|(a,b) and 
are not known analytically; nonetheless their functional form is similar to that 



of the Kolmogorov flow (5.5) 



To understand the diffusive 0{e^''^) scaling in the cigars and to determine 
something of the local structure of the fine-scaled A field the following strat- 



egy is adopted: solve the equations (6.10) and (6.11) by integrating along 



characteristics of A given locally by (6.18) and reconstruct A via (6.13). As 



the characteristics of A approach the origin and are squeezed along the out- 
going separatrix, given hj ( = 0(1), p = 0, high gradients build up and the 
terms in e that were earlier neglected increase: when these come into balance 
with the terms we have retained, we reach the scale at which diffusive effects 
become important, fixing the width of the dissipative regions. 

There are two problems with this approach: first that the incoming values of 
a and c are determined by the outer solution and links to other cigars: as 
this is beyond what can be addressed analytically, unknown functions have to 



be introduced. Secondly, even with the simplified, general local form (6.18) 



analytical calculations rapidly become unwieldy. The first problem will remain 
with us, but to ameliorate the second problem we simplify further and consider 
only the motion near to the origin, in which will simply take the field A to 
be, exactly, 

\ = {-ap,up,2aC), (6.19) 

in the local cylindrical polar coordinate system {p, </),()■ Here a and u are 
taken as constants, which we may identify as 



a = s'(0), u = n{0). 



(6.20) 



We also note from (6.19) that 

Vx A = (0,0,2w), A'^ = {uj^ + a'^)p'^ + Aa'^C =Aa'^C + 0{p'^). (6.21) 



Our strategy now is to solve the outer, diffusionless equations (6.10-6.11) for 
transport of a and c for the simplified form (6.19) of A. This is done exactly. 
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but then to see how large the neglected, diffusion terms are, we approximate 
by taking ( = 0(1) but p <^ 1, so our results are valid in the outer region, 
near to the origin, on the outward-going separatrix, as depicted schematically 
in figure [SJ Of course, by the time ( = 0(1) we are, strictly speaking, away 



from the stagnation point at the origin and the form (6.19) that we are using 
no longer applies. However the above form is sufficient to obtain the overall 
structure of the outer solution as the separatrix is approached, together with 
the scaling of the diffusive layer width. 



Equation (6.10) becomes 

A . Va = -AcoaC, (6.22) 

and letting t be a time parameter along characteristics, we integrate this in 
the standard way, with 

p = —crp, 9 = u, ( = 2a(, a = —Aua^, (6.23) 

and the solution in terms of initial conditions on a characteristic, 

p = Poe-"*, e = eo + ujt, C = Coe'"*, a = ao + 2u;Co(l-e2"*)- (6-24) 

If we suppose that we specify the incoming values of a on a surface p = po > 
(see figure [s]) with 

ao = a(po,^o,Co) = ^(^o,Co) (6.25) 

at t = 0, then we have the solution: 

a{p, e, = A[e + a-'u log(p/po), CpVpo] + 2a;C(pVPo - !)• (6-26) 

Here A gives the form of the field being carried in from the outer region, and 
we do not know much about it, except that it has 3-fold rotational symmetry 
(see figure |3] and figure 13 of CG). It is perhaps helpful to think of A as being 
some function of order unity with the appropriate symmetry, for example 
A = y4o + A3COS36'. 



Given a we can now reconstruct the appropriate part of A in (6.13). We have 



Va = {ua-^p-^Ae + 2CpPo\A^ + 2u),p-^Ae, p'po '^C + 2w(p'po ' " 1)) 
= a-^p-\uAe, aAg, -2uap) + 0{p), (6.27) 

and so 

Xa = A^^A X Va = 2A-^{-a, tu, 0)Cp'^Ae + 0{p^), (6.28) 

as p — )■ 0, where Aq denotes the derivative of A with respect to its first 
argument. Here we have obtained a component growing as p~^ which arises 
because of the incoming values of A on different characteristics being squeezed 
together. 
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6. 3 The effect of diffusive terms 



With this component of A in hand, we can now revisit the diffusive equation 



(6.7). We calculate 



V X A, = 2A-^(uj^ 



+ a')iO,0,a-')Cp-'Age + Oip 



,-l^ 



(6.29) 



and 



A ■ V X A, = AA-\u^ + a^)Cp-^Aee + 0(p-^). (6.30) 

This now has a p~^ growth, by virtue of differentiating A again. In equation 



(6.7) it is clear that the final eX ■ V x A^ term with diffusion will be the same 
order as the term A ■ V x A = 4:UjaC, we originally included, when ep"^ = 0(1). 
This gives the p = 0{^/e) scaling of the diffusive cigar width. Similarly at these 



values of p, in (6.8) and (6.9) the terms eV x \a become of similar magnitude 



to V X A (though note here that the Va terms are actually larger in magnitude 
at this point). 

This is the main part of the argument: although we have simplified by focusing 



solely on Aa in (6.13), consideration of the scalar c and component Ac does 



not affect the discussion, nor does Aa, given straightforwardly by 
Aa = A^^A X (V X A) = 2K-^uj{uj, a, 0)p, 



(6.31) 



and so is negligible. To check this we now look at the e = equation (6.11) 
for c and the corresponding component Ac. After a straightforward calculation 



(6.11) becomes 



A • Vc = -12A-V(a;2 ^ a^)^Ae + 0{p). 



(6.32) 



The key point is that the right-hand side is of order unity as p — )■ 0, as was 



the case for a in (6.22). Thus without solving the equation in detail, it is clear 



that the solution analogous to (6.26) for a will take the form 



c{p, e, = c[e + ioa-' iog(p/po), CpVpo] + Cpi(p, e, c), 



(6.33) 



where C{9o, Co) gives the incoming values of c on the surface p = po, as before 
and the particular integral Cpj involves A but is of order unity as p — )■ 0. 



Now when we reconstruct A via (6.13), the component cA 



0(1) along 
streamlines will be subdominant to the component A x Va = 0{p~^), the 
inverse power of p arising from taking the gradient. Thus our focus on Aq in 
the above discussion of the diffusive breakdown of the outer solution is justified 
and we have 



A = A, + 0(1) = 2A-^(-a, oo, 0)Cp-M, + 0(1) 



(6.34) 



16 



as p — )■ on the outgoing separatrix. As a by-product of our calculations we 
observe that the small-scale field A will show components Aq perpendicular to 
streamlines that diverge as p~^ as the separatrix is approached from the outer 
solution. These will peak at levels A = 0(e~^/^) when diffusive suppression 
begins to occur at scales p = ^. This is in keeping with the scalings seen 
by CG, who note that A_ = e\ = 0[^\fe) near the separatrices (their section 
3.2.1, figures 13 and 16). In view of the cos 36' dependence of the leading 
field identified here, this component must go to zero on the axis itself and is 
presumably strongly suppressed by diffusion. Thus we cannot make a detailed 
link with figure |4j the field here originates with the mean component of A, 
independent of ^, for which the onset of diffusion will be delayed until smaller 
values of p. This also presumably explains the structure seen in figure |3] (most 
clearly in (b)) or figure 13 of CG, with three incoming sheets of field merging 
in an axisymmetric 'collar' at smaller values of p. In this way, there could be 
several nested boundary layers along the separatrices in the limit £ — ?■ 0. 



7 Existence of weak steady-state solutions 



We consider now the system of equations (6.2) and (6.3) in Elsasser variables, 
together with the solenoidality conditions 

V-A = V-A = 0. (7.1) 

In this section we define weak solutions to these equations and formally prove 



their existence, adapting the approach of Ladyzhenskaya (1969). In the next 
two sections we will show that the weak solutions are classical smooth func- 
tions, satisfying the equations at any point in space. 

We start by recalling some definitions. Consider the class of functions whose 
domain is the periodicity cell T^ = [0, 27r]'^. The norm in the Lebesgue space 
Lp(T^) is defined, for p > 1, as 

Since in the above-mentioned class X — V^ is a positively defined self-adjoint 
operator (where X is the identity), whose eigenfunctions are Fourier harmonics, 
we can define in the usual way the powers (X — V^)" for an arbitrary real a, 
by considering Fourier series. For r ^ B? and any 

(7.3) 
(X - V')"$ = ^ (1 + H'T^r^ e'"-". (7.4) 
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The Sobolev space Wp{T^) is defined for p > 1, as tlie closure in tlie norm 



l$l 



s,p 



fX- V^)"/^$| 



(7.5) 



of the set of infinitely smooth periodic functions, whose domain is T^. (Evi- 
dently, Lp(T^) = Wp{T^).) We will work in the subspace of zero-mean vector 
fields, in which the operator — V^ can be used instead of X— V^ in these defi- 
nitions. In particular, we define (without introducing a new notation) a norm. 



equivalent to (7.5), in the subspace of zero-mean fields in W^(T^) as 



l$l 



s,p 



.V^)''/^$| 



(7.6) 



Since the Laplacian is a self-adjoint operator, in the important particular case 
p = 2 this implies 



l$l 



s,2 



^■(-V'Y^dV. 



f3 



(7.7) 



We will employ the following: 



Embedding theorem (see Bergh & Lofstrom (1976), Taylor (1981) and ref- 
erences therein). 

(i) For s > N/p, W^{T^) C ^(T^). 

(ii) For < s < N/p and q = Np/{N - ps), W.^{T^) C ^^(T^) (in particular, 

||$|| < c ll$ll ) 

We will show in the remainder of this section that for any space-periodic 
forcing / from the Lebesgue space L2(T^), the system of equations (6.2), (6.3) 



and (7.1) has at least one weak space-periodic solution from the Sobolev space 
14^2^ (T^). The assumption that the box of periodicity is the cube T'^ = [0, 27r]^ 
is technical: our arguments can be repeated almost literally for the case of 
an arbitrary parallelepiped of periodicity. Note that in this and the following 



sections we do not restrict ourselves to the Kolmogorov forcing (2.4); higher 
regularity of / will be required in ^ 

Consider then, the set of infinitely smooth solenoidal zero-mean periodic func- 
tions, whose domain is the periodicity cell T'^ = [0, 2tt]^, and denote by H its 
closure in the Sobolev space W2{T^). A pair of vector fields A ETi, A G H is 



a weak solution to the system (6.2), (6.3) and (7.1), if the integral identities 
and 



1-3 



'^i:l^-l^.((A-v)A-/).* 



fc=i 



dxk dxi 



dV = 



(7.9) 



hold true for any vector field $ G "H. (If A and A are smooth, these identi- 
ties immediately follow from (6.2) and ( |6.3 ).) By Holder's inequality and the 
embedding theorem, for any function / G VF2^(T^), 



^<||^||l/4„,„3/4 



<Cl||/llfll/llS'<^l||i 111,2, 



(7.10) 



where Ci is a constant independent of /. Consequently, the Cauchy-Bunya- 
kowsky-Schwarz inequality implies that the integrals involving nonlinear terms 
admit the bounds 



/ {{X-V)\)-^dV = f Y.^kA-j—dV 
Jj3 Jri fH dxk 



k=l 
3 3 

- 5UI ll^felUIIAjlU 
j=i k=i 



9$. 



dxk 



(7.11) 

^ C*2||A||i^2||A||i^2||^||l,2; 



C2 being a constant independent of A, A and $, and similarly 



/ ((A-V)\)-^dV 
Thus the integrals are well-defined. 



< C'2||A||i,2||A||i_2||^| 



1,2- 



(7.12) 



Consider the scalar product in Ti 

[$1,$2] = 



35$! 5$2 „. 



(7.13) 



Integrating by parts we recast the identities (7.8) and (7.9) in an alternative 



form involving the scalar product (7.13): 



and 



Here 



[A-^(A,A)-/,$]=0, 



[eX-A{A,X)-f,^] = 0. 



/ = -(v^)-^/ 



(V ) denoting, as usual, the inverse Laplacian, and 

^(A,A) = (V2)-V((A-V)A) 



(7.14) 

(7.15) 
(7.16) 

(7.17) 



is a bilinear operator, where V is the projection onto the subspace of solenoidal 
vector fields. (In fact, for the Kolmogorov forcing / = /, but in what follows 
we do not employ this equality.) 
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Using (7.7) for s = 1, we find 



\\A{\, A) II?,, = - / V{{\ ■ V) A) . {V'r'ViiX ■ V) A) dV 



^ — 1 ; 1 -^ -II- 



d^ 



dxjdxk 



(A,A) 



dV. 



For any 



{VT V 



d' 



dxjdx 



$ 



jUXk 






n^O 



I 19 11 19 

nr / \n\ 



and therefore 



(V^)-ip 



a^ 



dxjdxk 



$ 



< ll$l 



(7.18) 

(7.19) 
(7.20) 

(7.21) 



Now we develop (7.18), using Holder's inequality and the embedding theorem, 



3 3 3 

\\A{\ X)\\l^ <T.mi IIAfclUIIA/IUIIA.AIh < C3||A||?,2||A||?,2, (7-22) 
j=i fc=i 1=1 

which shows that A : Ti^'H ^ Ti. 

Thus, we have shown that for A G H and A G "H the first factors in the 



scalar products in the right-hand sides of (7.14) and (7.15) belong to "H. Since 



smooth vector fields are dense in "H in the norm induced by the scalar product 



(7.14) and (7.15) are equivalent to equations 



A-^(A,A)-/ = 



and 



\-e'\A{\A) + f) = 0, 
respectively, understood as equalities in Ti. 



(7.23) 
(7.24) 



The existence of solutions to the system (7.23), (7.24) is guaranteed by the 



Leray-Schauder principle (see Leray & Schauder (1934) and Ladyzhenskaya 



(1969)) under two conditions: 



(i) The operator B : U ^U ^ U ^U defined as 

S(A,A) = (^(A,A),^(A,A)/£) 



(7.25) 



is compact, i.e. i3(A„, A„) is a strongly converging sequence in "H ® "H for any 
sequence (A„, A„) weakly converging in "H ® "H. 
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(ii) Any solution to the set of equations 

A-/i(^(A,A) + /) = 0, A-/i£-i(^(A,A) + /) = (7.26) 

belongs to a ball in H "H of a radius independent of yU for < yU < 1 . 

The proof of (i) relies on the embedding theorem for Sobolev spaces, whereby 
the embedding W^{T^) -^ Lg{T^) is compact for g < 6, i.e., ||A" - A"||g + 
II A" — A^llg — 7- for m,n —> oo, for any sequence (A", A") weakly converging 
in 14^2^ (T^) 14^2 (T^). It is enough to prove that ^(A", A") converges strongly 
in 7i. For any $ G H, 



[^(A",A")-^(A'",A'"),$] = / ^(A^A"-A^A™)-^— rfV^ 



3 
fc=l 



9$ 



ltKi^"-^n-^/v^lhx 



A^)A" 



dxk 
dxk 



dV. 



(7.27) 
(7.28) 



Hence, by the same arguments as were used to derive (7.11), we obtain 

|[^(A",A")-^(A'",A'"),$]||i,2<C4(||A"||i,2||A"-A™||4+||A"-A"^||4||A"^||i,2)||$||i,2. 

(7.29) 
Here $ G "H is arbitrary; letting $ = ^(A", A") - ^(A™,A'"), from this 
inequality we deduce 



M(A",A")-^(A™,A^ 



11,2 



<a(IIA"-A" 



+ ||A"-A'"||4), (7.30) 



where the constant C5 is independent of m and n, since weak convergence of 
(A", A") in W^{T^) ® W^{T^) implies the uniform boundedness of || A"|| 1,2 and 
||A"||i,2. Thus we have established that ||^(A", A") - ^(A™, A™)||i,2 -^ for 
m, n — )■ 00, as desired. 



To prove (ii), we consider the problem (7.26) in the form of integral equations. 



analogous to (7.8) and (7.9), 



T3 



3 P)\ f)^ \ 



(7,31) 



which are satisfied for any $ G "H. Let $ = A in (7.31 ) and $ = A in ( |7.32[ ). 
Due to solenoidality of A and A the nonlinear terms vanish, and we find from 



the identities (7.31) and (7.32) 



[A,A]=/i/ /■Arfy</iC6||A||2||/||2 => ||A||i,2</iC7||/||2, (7.33) 
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[A, A] = fie-' [ fXdV< iiE-'Ce\\Xh\\f\ 



=> ||A||i,2 <^£ ^CrWfy, 

(7.34) 

since the norm, induced by the scalar product [■, ■] in T-L, is equivalent to the 
norm (7.7) in 14^2^ (T^). These inequalities establish the existence of a weak 



solution A E Ti, X E Ti to the problem (6.2), (6.3) and (7.1), admitting the 
bounds 

(7.35) 



lAII <^ n W -fW ll\ll <^ n c^^W -FW 

K»-||l,2 S tv7||J||2, ll'^l|l,2 S ^7^ \\J\\2- 



Here the constant Cy is absolute: it is independent of the solution A and A, 
the forcing /, and the parameter e. 



8 Bounds for weak solutions in W^(T'^) 



In this section we obtain bounds for the norms of the weak solution, whose 
existence we have established in the previous section, in the Sobolev spaces 
W2 (T^) and iy|(T'^). The forcing / is assumed here to belong to the Lebesgue 
space L2(T^), as in the previous section. 



Consider the Fourier series 



A 






U) 



and smooth vector fields 



AI 



n=/=0, \n\<M 



en. 



1.2) 



Scalar multiplying in L2(T^) (7.23) by (—V^)^/^ A G "H, using self-adjointness 
of the Laplacian, solenoidality of A and hence of A , and orthogonality of 
potential and solenoidal fields in L2(T^), we obtain 



/ A ■ (-V')^/^A*^ rfV + / (A-V)A-(-V')^/^A*-'rfr 

Jf'i 



Note that 



-V^)i/^((-V^)i/^A^0ll2 = ||A'"||i,2<||A||i,2, 



and hence {-V^f/^X^^ e Wl''^{T^), and by part (ii) of the Theorem 

||(-V2)V4A^^||3<Ci/2,2||A||i,2. 



U) 



U) 



^.5) 
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This, together with (7.35), imphes a bound for the first integral 



/ (A-V)A-(-V2)i/^A^^rfF 

JT3 



< ||A||J|VA||2ll(-V^)^/^A^^| 



Also, 



<Ci/2,2W\l2\\\\\i,2<Cs/e. (8.6) 

/ f-{-y'f"K'UV <||A||i,2||(-V2)^/^/||2<C9, (8.7) 



and hence we find from the identity (8.3) 



||(_y2)5/8^||2 = sup / A ■ {-V^f'^X^dV < Cio^^^ (8.8) 



Similarly, (7.24) yields 



-V'f/'XWl < Cios-'. (8.9) 



(The constant Cio in (8.8) and ( 8.9[ ) is independent of £ < 1, but depends on 
the norm ||/||2 of the forcing / G L2(T^).) We have therefore demonstrated 
that A G W2^\T^) and A G W^^^T^). 

Consequently, {-V^/^A^' G W^^^T"") and {-V^^^A'' G W^^^T""). Using 
part (ii) of the theorem, we find 

||9A/axfc||i2/5 < Cy,,2\\{-^'y^'{dA/dxk)h < Ci/,,2\\{-V'f'Ah (8.10) 
and 

||A||i2<C5/4,2||(-V2)5/«A||2. (8.11) 

Similarly, 

\\dX/dXk\\l2/, < Ci/,,2\\{-V'f'X\\2, l|Al|i2 < C5/4,2ll(-V2)'/'A||2. (8.12) 



Scalar muhiplying in L2iT^) (7.23) by {-WyA'', we obtain 



[-V')A'' \\i + (A ■ V)A ■ (-V')A^-' dV= / ■ (-V')'A^^^ dV. (8.13) 

JT3 JT3 



By Holder's inequality. 



/ (A-V)A-(-V')A*'rfl^ 

JT3 



< IIAIUIIVA 



1 12/5 I 



-V2)A*^||2 



(8.14) 



and hence from (8.13) and (8.10-8.12) 



-V^)A^||^ < C^^\\{-Vr^'X\m-V'f/'A\\l 

+ i||(-V2)A^^||^+||/||2||(-V2)A| 



2, 



^.15) 
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whereby 



-V^)A||^ = supl|(-V^)A^'^||^<C7l2£-^ 

A/ 



(8.16) 



The same operations apphed to (7.24) yield 



-V')A||2 < Ci2e-\ 



(8.17) 



Thus we have demonstrated that A and A belong to Wg ('^ )• The constant 



Ci2 in (8.16) and (8.17) is independent of the solution A, A and the small 



parameter e, but depends on the norm ||/||2 of the forcing / G L2(T^). 



9 Smoothness of weak solutions 



Steady-state hydrodynamic and MHD problems are drastically different from 
the evolutionary ones in that one can incrementally establish the smoothness 
of their solutions together with the derivatives of arbitrary order (provided 



the forcing / is sufficiently smooth). In this section we use (7.23) and (7.24) 
to show by induction that A and A, whose existence we have ascertained in 
^\ are in fact smooth vector fields and therefore constitute a classical space- 



periodic solution to equations (6.2), (6.3) and (7.1) 



We assume now that A G W2^{^^) and A G W2'^{T^) (which is equivalent to 

||(-V2)'=A||2 + ||(-V2)'=A||2 < oo) and / G W^^iJ^) for A; > 1, and show that 
A G W^^^^it^) and A G W^'^+^iT'^). 



Scalar multiplying in L2(T^) (7.23) by (— V^)^'^"'"^A*^, using self-adjointness 
of the Laplacian, solenoidality of A and orthogonality of potential and 
solenoidal fields in L2(T^), we obtain 



Jt3 ^"'t3 \^ dxj J 






JT3 



(9.1) 



and thus 






-V')''"'7;f-(A-V)A 

dXj 



(9.2) 



+ ||(_V2)'=-V2j|| ^ ||(_v2)fc+l/2^M| 
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implying 



|A||2fe+l,2 = sup 
M 



._y2)/c+l/2^Af||^<^ 



-V 



2\fc-l 



d_ 
dxj 



:a-v)a 



+II/I 



2fc-l,2- 



(9.3) 



We need therefore to check that the norms in the sum at the right-hand side 
of this inequahty are bounded. By part (i) of the Theorem, the assumption 
A G W2''(T^) and A G 1^2^'' (T^) imphes that A and A and their derivatives 
of order up to 2A; — 2 are continuous (and hence uniformly bounded) vector 
fields in T'^. By the standard formula for derivatives of products, 



72\fc-l 



dXr 



(A-V)A 



(9.4) 



is a linear combination of products of derivatives 



9^U„ 



d^'A 



9"iXi9"'2X2(9"3X3 9'^iXi9"2X2(9"'3X3 



(9.5) 



where Ni = n\ + ni + ni, < Ni < 2k - 1, N2 = 2k - Ni. Thus, each 
of the terms is continuous and bounded, except maybe those for Ni = 0, 1 
or 2k — 1. If A^^i = 0, or A^^i = 1 or 2/;; — 1 for k > 1, one of the factors 
is continuous and the second one is known to belong to L2(T^), and hence 



their contributions to the right-hand side of (9.3) are finite. The remaining 
possibility is A; = A^^i = A'^2 = 1, but in this case both factors belong to 
L4(T=^) because A G Wi{T^) and A G W^{T^) by the results of the previous 
subsection, and again the respective norms are bounded. 

We h ave thus established that A G W2'^'^(T^). By similar arguments from 
(fl^ we find A G Wi''+\T^). 



To proceed, we scalar multiply in L2iT^) (7.23) by (-^^y^+^A^, and obtain 



/ A ■ (-V')''+'A^^ dV + / ((-V')'=(A ■ V)A) ■ (-V2)'^+^A^^ dV 

= I i-V^ff-i-V^f^^A'^dV, (9.6) 

JT3 



and therefore 



lAI 



2fc+2,2 



sup 

M 



2\k, 



_^2).+i^.w||^ < ||(-V^)'^(A ■ V)A||2 + 11/ 



2fc,2- 



(9.7) 



Since A G Wi^^\T^) and A G Wi''^\T^), by part (i) of the Theorem any 
derivative of A and A of order up to 2A; — 1 is continuous in T^. Hence, in the 



expansion of (A ■ V)A in a linear combination of products (9.5) each term is 
either continuous, or a product of a continuous function by a function from 
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L2(T^). Thus, (19/7|) demonstrates that A G Wf^\T-^). Similarly ( |7.24[ ) yields 
A G W2''~^'^ (T^) , concluding the demonstration. 

Thus mathematical analysis of the problem yields both good and bad news. 



The good news is that the problem (6.2), (6.3) and (7.1) necessarily has at 



least one classical solution, meaning in our case infinitely differentiable at each 
point. Nothing is known about the number of solutions except that it is strictly 
positive, nor is stability of any of the MHD steady states guaranteed. The bad 
news is that the bounds for the solutions and their derivatives rapidly degrade 
as e — > 0. In particular, the inequalities that we have derived are insufficient 



to claim that eV^A — )■ 0: the relevant L2 bound we have derived is (8.17) 



It is interesting to compare this general result, that is for a general forcing, 
with our nu meric al study of dissipative regions in the Archontis case. For 
example, in ^6.3 we find peak values A = 0{£^^^'^) on scales of order e^/^ 
indicating the scaling eV^A = 0{e^^^'^). Note that these anomalously high 
values are concentrated only in cylindrical cigars about the separatrices, of 
radius 0(e^/^) and so occupy an 0(e) volume of space. This results in the 
estimate ||A||2,2 = 0{e~^), which is a significantly milder singularity than the 



one suggested by our bound (8.17). The high values have a negligible impact 



on the energy spectrum, there being no peak visible at small scales. This can 
probably explain the gap between the 'worst case scenario' predicted by the 
rigorous mathematical analysis of the problem and the numerical results: the 
Sobolev norms, that we have used, prove inefficient in controlling formation 
of singularities in localised regions, because they are of inherently integral 
nature. We should also note that our simulations would not be able to resolve 
structure on scales much smaller than 0(e^/^). 

The apparent deterioration of the derivatives of the solution with their order 
can be a spurious artefact due to imperfection of the proof (which is especially 
possible in view of the generality of our arguments — at no point in §^|3-[9]we 
have made use of the fact that the dynamo that we consider is powered by the 



Kolmogorov forcing (2.4)), or a real attribute of the solutions. It is likely that 



for some forcing in (6.2) and (6.3) the worst case scenario suggested by these 



bounds is indeed realised: they are based on the norm bounds provided by 
the embedding theorem, which are sharp. In any case, this indicates that any 
naive approach to the study of the limit £ — )■ (for a general forcing) whereby 



the diffusive term in (6.3) is just discarded, is likely to be erroneous; this can 



only be done in the region outside dissipative structures. In the absence of 
the dominant elliptic operator, the equations obtained in this way are not in 
general guaranteed to have solutions. When they exist, the solutions are likely 
to develop singularities at some points or on certain manifolds, or possibly on 
sets of a more complex structure. Note that locally the existence of solutions 
is not a problem: the difficulty is in gluing together patches of such solutions. 
The fast dynamo problem embodies a similar difficulty, with small scales of 
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magnetic field occurring in 0(1) volumes of space, though with the fields 
concentrating on multifractal sets ([Childress fc Gilbert 1995). Note however, 



that the L2(T^) norm of A is uniformly (over e) bounded, so the singularities 
are likely to be more pronounced in the derivatives of the solution, rather than 
in the solution itself. 



10 Discussion 



We have presented investigations into the structure of the magnetic field and 
flow in the equilibrated regime of the Archontis dynamo. Because of the highly 
three-dimensional nature of the system, application of the available analytical 
tools yields only rough results of limited value, and we lack any kind of com- 
plete solution. What we have done is first to extend the range of diffusivities 
e over which the saturation mechanism operates to give the steady state with 
nearly aligned fields. We have also classified the symmetries of these flows 
and measured the field structure on the separatrices, home of the cigar-like 
dissipative regions. 

Then, using basic analytical tools, we have investigated the scaling of diffusive 
terms near the separatrices. Here at leading order the field A = £:~^ A_ that en- 
ters from the body of the flow is transported along characteristics of A = A+. 
Where these characteristics come together in the compressive flow at the stag- 
nation points, where trajectories spiral in, large gradients in A are generated, 
and diffusive terms enter the problem on scales of e^/^ as found by CG. In more 
general flows we may expect a similar behaviour, with regions of heightened 
dissipation localised at points where A = and along the unstable manifolds 
of such points. Of course in the Archontis example the unstable manifolds 
link the stagnation points and so the topology here is very simple and the 
dissipative regions very small, of order 0{6) in volume: in other cases they 
may wander through the three-dimensional space, giving a picture of much 
greater complexity, as could be occurring in examples in ICameron & Gal- 



loway (2006b). Again wider regions of dissipation, perhaps dense in the space, 
could occur if examples exist where A has no stagnation points; unfortunately 
the form of A is not under our control except where strongly constrained by 
symmetries. In order to cope with unknown levels of geometrical complexity, 
an approach based on functional analysis is appropriate, and this is the final 
part of the paper, in which the existence and smoothness properties of steady 
solutions are established. 



An analogy of the naively truncated equations (namely (6.2, 6.3) with e = 0) 
with the Euler equation, which is the subject of intense research, is instructive. 
A method for the investigation of the evolutionary Euler and Navier-Stokes 
equations consists of the introduction into the equations of new regularising 
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terms, such as e{—V'^)"u or e{—V'^Y{du/dt). It has been known for decades 
that for a > 5/4 solutions to the regularised equations are infinitely differen- 
tiable at any t > 0; for /3 > 1/2 and (3 > 5/6 one can prove analyticity, at 
any t > 0, of solutions to the Navier-Stokes and Euler equations, respectively 



( Zheligovsky , 2010). For any a or (3 below the respective thresholds, the prob- 
lem is as difficult as the one for the original equation. When the limit £ — > is 
considered, the results are so far inconclusive. One can only show that there 
exist sequences e^ — >• such that solutions for these Ek converge to a weak 
solution to the non-regularised equation, and either the limit weak solution 
is unique for all such sequences, or there exists a continuum of weak solu- 
tions. Whether for e — )• singularities develop in derivatives of the regularised 
solutions, and how strong they are if they develop, remains unknown. 

The difficulties arise in the general theory, because the bounds for solutions are 
singular in e as e — )■ 0. Here the analogy with the Archontis dynamo problem 
crystallises: in the Archontis problem the diffusive terms can be regarded as a 
regularisation of the naively truncated diffusionless problem, and we need to 
find out what happens when the regularisation parameter e tends to zero. (In 
the diffusionless, i.e. non-regularised, case it is unclear whether weak steady 
solutions exist.) We note that the analogy may work both ways: the asymptotic 
analysis near the separatrix in the Archontis dynamo (which we present in §[|5] 
and [6]) may contain clues to what happens in solutions to the regularised Euler 
(or even Navier-Stokes) equations in the limit e — )■ 0. Unfortunately, the clues 
are well hidden, because the regularising term in the Archontis problem is 
of a different structure, and a very specific symmetric steady solution to the 
general system of MHD equations is considered. 



Besides further attempts to carry out an asymptotic analysis of equations (6.2) 



and (6.3) and their evolutionary versions, a number of other directions could be 
pursued in the future, for example investigating time-dependent modifications 
to the steady Kolmogorov forcing used here, or studying the evolution of 
superposed large-scale fields and corresponding non-helical transport effects. 



as in the recent work of Sur & Brandenburg (2009). 
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